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We describe and analyze an efficient register-based hybrid quantum computation scheme. Our 
scheme is based on probabilistic, heralded optical connection among local five-qubit quantum regis- 
ters. We assume high fidelity local unitary operations within each register, but the error probability 
for initialization, measurement, and entanglement generation can be very high (~ 5%). We demon- 
strate that with a reasonable time overhead our scheme can achieve deterministic non-local coupling 
gates between arbitrary two registers with very high fidelity, limited only by the imperfections from 
the local unitary operation. We estimate the clock cycle and the effective error probability for 
implementation of quantum registers with ion-traps or nitrogen-vacancy (NV) centers. Our new 
scheme capitalizes on a new efficient two-level pumping scheme that in principle can create Bell 
pairs with arbitrarily high fidelity. We introduce a Markov chain model to study the stochastic 
process of entanglement pumping and map it to a deterministic process. Finally we discuss require- 
ments for achieving fault- tolerant operation with our register-based hybrid scheme, and also present 
an alternative approach to fault-tolerant preparation of GHZ states. 



I. INTRODUCTION 

The key challenge in experimental quantum informa- 
tion science is to identify isolated quantum mechanical 
systems with good coherence properties that can be ma- 
nipulated and coupled together in a scalable fashion. Re- 
cently, considerable advances have been made towards 
interfacing of individual qubits in the optical and mi- 
crowave regimes. These include advances in cavity QED 
[U [2] as well as in probabilistic techniques for entangling 
remote qubits [H [H 0. At the same time, substantial 
progress has been made towards the physical implemen- 
tation of few-qubit quantum registers using systems of 
coupled trapped ions [H |71 [8] , neutral atoms [9^ , or solid- 
state qubits based on either electronic and nuclear spins 
in semiconductors [lOl [TTJ [12] or superconducting islands 
[T3l[14]. 

While the precise manipulation of large, multi-qubit 
systems still remains an outstanding challenge, various 
approaches for connecting such few qubit sub-systems 
into large scale circuits have been investigated |5l [151 
[161 113 [HI- These studies suggest that hybrid schemes, 
which benefits from short range interactions for local cou- 
pling and (optical) long range interactions for non-local 
coupling, might be an effective way toward large scale 
quantum computation: small local few-qubit quantum 
systems may be controlled with very high precession us- 
ing optimal control techniques [I9l [20] , and in practice it 
may be more feasible to operate several such small sub- 
systems compared to the daunting task of high-precession 
control of a single large quantum system with thousands 
of qubits. Optical techniques for quantum communica- 
tion can then be used to connect any two sub-systems. 
For example, we may directly transfer a quantum state 
from one sub-system to another via an optical channel 
[2T] (see Fig. [l^), which immediately provides an effi- 
cient way to scale-up the total number of the physical 



qubits we can manipulate coherently. In particular the 
use of optical means for connecting different subsystems 
has the advantage that it allows for fast non-local opera- 
tions over large distances. This is advantageous for quan- 
tum error correction since the existence of such non-local 
coupling operations alleviates the threshold requirement 
for fault-tolerant quantum computation [22] , 

In practice, however, it is very difficult to have a perfect 
optical connection. In particular, there is excitation loss 
associated with the optical channel, due to scattering or 
absorption. For a lossy channel, it is therefore more desir- 
able to use it to generate entanglement between different 
sub-systems (see Fig. [T]3), rather than for direct state 
transfer. The entanglement generation is then heralded 
by the click patterns from the photon detectors. Such 
detection-based scheme is intrinsically robust against ex- 
citation loss in the channel, since it only reduces the 
success probability but does not affect the entanglement 
fidelity. This entanglement can then be used as a re- 
source to teleport quantum state from one sub-system 
to another [23 . More generally, entanglement provides 
a physical resource to implement non-local unitary cou- 
pling gates (such as CNOT gate) [241 [25l [H [27, . If there 
is only one physical qubit for each sub-system, so-called 
cluster states [28^ can be created based on probabilistic, 
heralded entanglement generation (see Ref. [5 and refer- 
ences therein). Such cluster states can be used for univer- 
sal quantum computation [29 . If there are two physical 
qubits available for each sub-system, cluster states can 
be created deterministically [30 ; meanwhile one can also 
use these two-qubit sub-systems to implement any quan- 
tum circuit directly [4]. 

Furthermore, realistic optical channel connecting sub- 
systems has other imperfections beside excitation loss, 
such as the distortion of the polarization or shape of the 
wave-packet. These imperfections will reduce the fidelity 
of the heralded entanglement generated between the sub- 
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FIG. 1: (Color online) Two schemes to couple different reg- 
isters, (a) Deterministic state transfer from one register to 
the other j21| . (b) Probabilistic entanglement generation, 
heralded by distinct detector click patterns [3] [31] [32]. In 
the text, we argue that probabilistic, heralded entanglement 
generation is sufficient for deterministic distributed quantum 
computation. 



systems. To overcome these imperfections in the channel, 
entanglement purification schemes have been proposed, 
which may create some high-fidelity Bell pairs from many 
low-fidelity ones 23 [34]. In particular, the entanglement 
pumping scheme originally presented in the context of 
quantum communication over long distances [SU [35l [36] 
provides a very efficient purification scheme in terms of 
local physical resources, and Diir and Briegel [15] first 
proposed to use such entanglement pumping scheme for 
quantum computation. In principle, the infidelity of the 
purified Bell pair shared by the sub-systems can be very 
low, and is only limited by the error probability from lo- 
cal operations. In Ref. [T5^ it was found that three auxil- 
iary qubits (requiring a total of five qubits including the 
storage and communication qubits) for each sub-system 
provide enough physical resources to obtain high fidelity 
entangled pairs via entanglement pumping. 

In order to implement the idea of distributed quantum 
computation using realistic optical channel and imper- 
fect operations, it is necessary to consider the following 
questions: What are the minimal local physical resources 
needed for robust entanglement generation? What is the 
time overhead associated with entanglement generation? 
Can we extend the robustness to other imperfections, 
such as errors from initialization and measurement? 

Motivated by these considerations, we study the prac- 
tical implementation of robust quantum registers for scal- 
able applications. In Ref. [37 we have proposed an 
entanglement purification scheme that only requires two 
auxiliary qubits for robust entanglement generation. We 
have found that the time overhead associated with entan- 
glement generation ranges from a factor of 10 to a few 
100, depending on the initial and targeting infidelities. 
We have also suggested to use one more auxiliary qubit 
to suppress errors from initialization and measurement. 
Thus, our hybrid scheme also requires only 5 (or fewer) 



qubits with local deterministic coupling, while providing 
additional improvements over the protocol of Ref. [15j: 
reduced measurement errors, higher fidelity, and more ef- 
ficient entanglement purification. In this paper, we will 
provide a detailed discussion on the register-based, hy- 
brid quantum computation scheme presented in Ref. [37] . 

The paper is organized as follows. In Sec. [TTj we will 
introduce the concept of quantum register and discuss 
two experimental implementations. In Sec. |III[ we will 
review the idea of universal quantum computation based 
on two-qubit quantum registers. In Sec. |IV[ we will 
specify the error models for imperfections and provide 
the basic ideas underlying our robust operations. In Sec. 
[v] we will describe the robust measurement /initialization 
scheme. In Sec. |VI[ we will present our bit-phase two- 
level entanglement pumping scheme. In Sec. |VII[ we 
will introduce the Markov chain model to quantitatively 
analyze the time overhead and residual infidelity associ- 
ated with the stochastic process of entanglement pump- 
ing, and discuss further improvement upon our two-level 
entanglement pumping scheme. In Sec. VIII[ we will map 
our stochastic, hybrid, and distributed quantum compu- 
tation scheme to a deterministic computation model that 
is characterized by two quantities (the clock cycle and ef- 
fective error probability), and estimate the practical val- 
ues for these quantities. We will also consider the con- 
straint set by the finite memory lifetime and determine 
the achievable performance of hybrid distributed quan- 
tum computation. Finally, in Sec. [iXj we will discuss us- 
ing our hybrid scheme for fault-tolerant quantum compu- 
tation with quantum error correcting codes, and provide 
a resource-efficient approach for fault-tolerant prepara- 
tion of the GHZ states. 



II. QUANTUM REGISTER AND 
EXPERIMENTAL IMPLEMENTATIONS 

We define a quantum register as a few-qubit device 
(see Fig. [2|i) that contains one communication qubit (c), 
with a photonic interface; one storage qubit (s), with 
very good coherence times; and several auxiliary qubits 

(ai,a2, ), used for purification and error correction. 

A critical requirement for a quantum register is high- 
fidelity unitary operations between the qubits within a 
register. 

The quantum registers considered here can be imple- 
mented in several physical systems, but in this paper we 
shall focus on two specific systems where these considera- 
tions are particularly relevant. First, ion traps have been 
used to demonstrate all essential elements of quantum 
registers. (1) The ion qubits may play the role of com- 
munication qubits: they can be initialized and measured 
efficiently using optical pumping and cycling transitions 
respectively, and they can also be prepared in highly 
entangled states with the polarization of single photons 
[39| [40] . Very recently, entanglement generation between 
ion qubits from two remote traps has been demonstrated 
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FIG. 2: (Color online) Distributed quantum computer, (a) 
Illustration of distributed quantum computer based on many 
quantum registers. Each register has five physical qubits, in- 
cluding one communication qubit (c), one storage qubit (s), 
and three auxiliary qubits (01,2,3)- Local operations for qubits 
within the same register have high fidelity. Entanglement be- 
tween non-local registers can be generated probabilistically 
1311 ^2 . Devices of optical micro-eletro- mechanical sys- 
tems (MEMS) [38 can efficiently route photons and couple 
arbitrary pair of registers. Detector array can simultaneously 
generate entanglement for many pairs of registers, (b) An 
ion trap coupled to a cavity also provides a promising can- 
didate for distributed quantum computation. A single ion is 
resonantly coupled to the cavity and serves as the commu- 
nication qubit; while the others can be storage or auxiliary 
qubits. (c) Nitrogen-vacancy (NV) defect center in photonic 
crystal micro-cavity. The inset shows the atomic structure 
of the NV center J^, which forms a quantum register. The 
electronic spin localized at the vacancy is optically accessi- 
ble (measurement/initialization) and can play the role of the 
communication qubit. The nuclear spins from proximal ^^C 
atoms constitute the storage and auxiliary qubits, which are 
coherently controlled via hyperfine interaction and rf pulses. 



j4T| [42] . This experiment directly demonstrates the non- 
local connection required for our hybrid approach. Since 
the photon collection and detection efficiency is not per- 
fect, the entanglement generation is a probabilistic pro- 
cess. However, the entanglement generation is also a her- 
alded process, because different click patterns from the 
detectors can be used to identify each successful entan- 
glement generation. As we will discuss Sec. |III[ such 
probabilistic, heralded entanglement generation process 
is already sufficient to implement deterministic non-local 
coupling gates. (2) The ion qubits can be good storage 
qubits as well. Coherence time of approximately 10 sec- 
onds has been demonstrated in ion traps [43 , 44 , which is 
-|^q6~7 i^ij^gg longer than the typical gate operation time 
which is at the order of jas. Since fault-tolerant quan- 



tum computation only requires the coherence time to be 
approximately 10^ times longer than the gate time, the 
very long coherence of the ions provides new opportu- 
nities, such as performing non-local coupling gates with 
some extra time overhead. (3) Coherent manipulation of 
few ions in the same trap has been demonstrated j45l |46] , 
allowing gates to be implemented among the qubits in 
the register. (4) High fidelity operations between the ion 
qubits within the ion trap has also been demonstrated 

A second promising candidate for implementing quan- 
tum registers is the nitrogen vacancy (NV) centers in di- 
amond. Each NV center can be regarded as an ion trap 
confined by the diamond crystal, which can be treated as 
a single register. The qubits for each NV register consists 
of one electronic spin associated with the defect and sev- 
eral nuclear spins associated with the proximal C-13 nu- 
clei. The electronic spin is optically active, so that it can 
be measured and initialized optically. With optical cavi- 
ties or diamond-based photonic crystal micro-cavities [47] 
one could enhance photon collection efficiency towards 
unity. Furthermore, the electron spin can be coherently 
manipulated by microwave pulses [10]. The electronic 
spin is thus suitable as the communication qubit. The 
nuclear spins are coupled to the electronic spin via hy- 
perfine interactions. One can either use these hyperfine 
interaction to directly rotate the nuclear spins ^48| |49j , 
or apply radio-frequency pulses to address individual nu- 
clear spins spectroscopically [50 . These nuclear spins 
have very long coherence times approaching seconds [12] , 
and can be good storage and auxiliary qubits. Further- 
more, the optical manipulation of the electronic spin can 
be well decoupled from the nuclear spins [51^. It can be 
inferred from the recent experiment [12] that the fidelity 
of local operations between electronic and nuclear spins 
is higher than 90%. While the fidelity is still low for 
the procedures considered here, we believe that it can be 
significantly improved (to higher than 0.999) by optimal 
control techniques [20] [50] , such as composite pulses [52] 
and numerically optimized GRAPE pulses [53]. 



III. UNIVERSAL QUANTUM COMPUTATION 
WITH TWO-QUBIT REGISTERS: 
FUNDAMENTALS 

We now consider universal quantum computation 
via the simplest possible two-qubit registers (H [30] . 
Each register has one qubit for communication and the 
other qubit for storage. We can use probabilistic ap- 
proaches from quantum communication ([31] and refer- 
ences therein) to generate entanglement between com- 
munication qubits from two arbitrary non-local regis- 
ters. The probabilistic entanglement generation creates a 
Bell pair conditioned on certain measurement outcomes, 
which are distinct from outcomes of unsuccessful entan- 
glement generation. If the entanglement generation fails, 
it can be re-attempted until success, with an exponen- 
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FIG. 3: (Color online) Quantum circuit for non-local CNOT 
gate between two registers and . The circuit starts 
with the Bell state ^2 = (|00) + 7^2 for the com- 

munication qubits and (the left blue box). After local 
unitary operations within each register (the middle orange 
boxes), qubits and are projectively measured in Z and 
X basis, respectively (the right green box). According to 
Eq. ([1]), up to some local unitary gates (the right green box), 
this circuit implements the non-local CNOT gate to qubits 
and from two quantum registers. 



tially decreasing probability of continued failure. 

When the communication qubits (c^ and c^) are pre- 
pared in the Bell state, we can immediately perform the 
non-local CNOT gate on the storage qubits {s^ and 5^) 
using gate-teleportation between registers and R'^. 
The gate-teleportation circuit in Fig. [3] implements (be- 
fore the conditional Pauli operations) the following map 



1^+)^, - «0-^^«0^^^CNOT,.,,. |(/>),,^,. , 

(1) 

where \^^)^i^^2 = (|00) + |11)) /V^, CNOT^,^- is the 
controlled-NOT (CNOT) gate with the ith qubit as con- 
trol and the jth qubit as target, and = 0, 1 is the 
measurement result for qubit i from the circuit in Fig. |3] 

By consuming one Bell pair, one can implement any 
non-local controlled-U gate between two storage qubits 
[27] , as shown in Fig.|4] Since operations on a single qubit 
can be performed within a register, the CNOT operation 
between different quantum registers is in principle suf- 
ficient for universal quantum computation [4 . Similar 
approaches are also known for deterministic generation 
of graph states [3^ — an essential resource for one-way 
quantum computation [29 . 

We emphasize that deterministic entanglement gener- 
ation is not required, which opens up a wide range of 
possibilities of entanglement generation. It is experimen- 
tally challenging to achieve deterministic quantum state 
transfer directly [21 , but we are able to achieve the same 
task by probabilistic entanglement scheme and two-qubit 
quantum registers [4 . 




FIG. 4: (Color online) General circuit for non-local 
Controlled-U gate [27], with the storage qubit S'^ as control 
and the storage qubit as target. One Bell pair is consumed 
for this operation. 



IV. ERRORS AND IMPERFECTIONS 

In practice, the qubit measurement, initialization, and 
entanglement generation can be noisy with error prob- 
abilities as high as ~ 5%, due to practical limitations 
such as imperfect cycling transitions, finite collection ef- 
ficiency, and poor interferometric stability. As a result, 
there will be a large error probability in non-local gate 
circuits. In contrast, local unitary operations may fail in- 
frequently {pl ^ 10~^) when quantum control techniques 
for small quantum system are utilized |19| i20j . We now 
show that the most important sources of imperfections, 
such as imperfect initialization and measurement errors 
for individual qubits in each quantum register, and en- 
tanglement generation errors between registers, can be 
corrected with a modest increase in register size. We de- 
termine that with just three additional auxiliary qubits 
and high-fidelity local unitary operations, all these er- 
rors can be efficiently suppressed by repeated quantum 
non-demolition (QND) measurement [54 and entangle- 
ment purification [35l |36] . This provides an extension of 
Ref. [15 that mostly focused on suppressing errors from 
entanglement generation. 

We will use the following error model for the entire 
paper: (1) The imperfect local two-qubit operation Uij 
is 

U.jpUl ^ (1 - pl) Ui^pU}^ + ^ TV,, [p] ® li,- (2) 

where Tr^j [p] is the partial trace over the qubits i and j, 
and lij is the identity operator for qubits i and j. This 
error model describes that with a probability I — Pl the 
gates perform the correct operation and with a probabil- 
ity Pl the gates produce a complete random output for 
the two involved qubits [71]. (2) The imperfect initial- 
ization of state |0) will prepare a mixed state 

po = (l-P/)|0)(0|+p/|l)(l|, (3) 

which has error probability p/, i.e., it prepares the wrong 
state with a probability pj. (3) The imperfect measure- 
ment of state |0) will correspond to the projection oper- 
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(4) 



This operator describes that a qubit prepared in state |0) 
or |1) will give rise to the opposite measurement output 
with the measurement error probability pm- (4) Finally, 
the entanglement fidelity for a non-ideal preparation of 
state |$+) is defined as 



F = 



(5) 



and the infidelity is just 1 — F. The fidelity does, how- 
ever, not completely characterize the produced entangled 
state. Depending on the exact method used to generate 
the entangled state, one can in some situations argue that 
the error will predominantly be, e.g., only a phase error 
|55] , whereas in other situations it will be a combi- 
nation of phase and bit fiip errors (see ^31^ and references 
therein). Below we shall therefore both consider the sit- 
uation where we only have a dephasing error as well the 
situation, where we have a more complicated depolariz- 
ing error (exact definition given later). As we shall see, 
the knowledge that the error is of a particular type (e.g. 
only dephasing error) provides a significant advantage for 
purification. 

We will also assume a separation of error probabilities: 
any internal, unitary operation within the register fails 
with extremely low probability, p^, while all operations 
connecting the communication qubit to the outside world 
(initialization, measurement, and entanglement genera- 
tion) fail with error probabilities that can be several or- 
ders of magnitude higher. 



PL < Pi.Pm, I- F. 



(6) 



In terms of these quantities the error probability in the 
non-local CNOT gate in Fig. [3] is 



PCNOT (1 - F) + 2pL + 2pM., 



(7) 



because we use one entangled state, two local operations, 
and two measurements. In the next two sections, we 
will show how to use robust operations to dramatically 
improve the fidelity for these non-local coupling gates. 

Robust measurement can be implemented by repeated 
QND measurement, i.e., a majority vote among the mea- 
surement outcomes (Fig. |5|, following a sequence of 
CNOT operations between the auxiliary /storage qubit 
and the communication qubit. This also allows ro- 
bust initialization by measurement. High-fidelity, ro- 
bust entanglement generation is achieved via entangle- 
ment pumping [151 EH (Figs. |6] and [t]), in which 
lower fidelity entanglement between the communication 
qubits is used to purify entanglement between the auxil- 
iary qubits, which can then be used for non-local CNOT 
operations. To make the most efficient use of physi- 
cal qubits, we introduce a new entanglement pumping 
scheme. In our bit-phase two-level entanglement pump- 
ing scheme, we first use unpurified Bell pairs to repeat- 
edly pump (purify) against bit-errors (Fig. ^), and then 
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FIG. 5: The robust measurement scheme based on repeated 
quantum non-demolition (QND) measurements and majority 
vote. Each QND measurement consists of initializing, cou- 
pling, and measuring the communication qubit. The QND 
measurements are repeated 2m + 1 times using the same com- 
munication qubit. 



use the bit-purified Bell pairs to repeatedly pump against 
phase-errors (Fig. [7]3). 

Entanglement pumping, like entanglement generation, 
is probabilistic; however, failures are detected. In com- 
putation, where each logical gate should be completed 
within the allocated time (clock cycle), failed entangle- 
ment pumping can lead to gate failure. To demonstrate 
the feasibility of our approach for quantum computa- 
tion, we will analyze the time required for robust initial- 
ization, measurement and entanglement generation, and 
show that the failure probability for these procedures can 
be made sufficiently small with reasonable time overhead. 



V. ROBUST MEASUREMENT & 
INITIALIZATION 

In this section, we will analyze the robust measurement 
scheme based on repeated QND measurement, discuss 
the recent experimental demonstration of robust mea- 
surement in the ion-trap system, and present two ap- 
proaches to robust initialization. 

The measurement circuit shown in Fig. [5] yields the 
correct result based on a majority vote from 2m + 1 con- 
secutive readouts (bit- verification) . Since the evolution 
of the system (CNOT gate) commutes with the measured 
observable {Z operator) of the auxiliary /storage qubit, it 
is a quantum non-demolition (QND) measurement [54J, 
which can be repeated many times. The error probability 
for such a majority vote measurement scheme is: 



/2m + l\ . ^ .J . 2m + l 

j=m-\-l 



PL (8) 



where the last term conservatively estimates the prob- 
ability for bit-fiip error of the auxiliary /storage qubit 
during the repeated QND measurement. For simplic- 
ity, we will use Eq. ([s]) for our conservative estimate of 
error probability for repeated QND measurement. Sup- 
pose Pi = Pm = 5%, we can achieve sm ^ 8 x 10~^ by 
choosing m* = 6 for = 10~^, or even sm ^ 12 x 10~^ 
by choosing m* = 10 for p^;, = 10~^. For convenience 
of discussion, we shall add £m to the set of imperfection 
parameters: {1 — F^pj^pm^Pl^^m)- The time for robust 
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measurement is 



tM = (2m + l)(ti + tL+tM) 



(9) 



where t/, t^,, and tM are times for initialization, local 
unitary gate, and measurement, respectively. 

Measurements with very high fidelity {em as low as 
6 X 10~^) have recently been demonstrated in the ion- 
trap system [56 , using similar ideas as above. There are 
several possibilities to further improve the performance 
of the repeated QND measurement. (1) We may use 
maximum likelihood estimate (MLE) to replace the ma- 
jority vote for repeated measurements with multi-value 
outcomes (e.g., fiuorescent intensity) [56^. (2) We may 
keep updating the error probability using MLE after each 
measurement. Once the estimated error probability is be- 
low some fixed error rate, we stop the repetition of the 
QND measurement to avoid errors from redundant oper- 
ations [56 . (3) We may use the implementation of CNOT 
gate that has small/ vanishing bit-fiip errors to the con- 
trol qubit, which will reduce/eliminate the last term in 
Eq. (||. 

The robust measurement scheme also allows to achieve 
robust initialization by measurement, i.e., by measuring 
the state of a qubit with the robust measurement scheme, 
we initialize the qubit into the result of the measurement 
outcome with an effective initialization error 



(10) 



Besides the above measurement-based scheme, we 
may achieve robust initialization using verification- 
based scheme [72]. For clarity, we shall assume the 
measurement-based initialization [Eq. ([lO|] for the rest 
of the paper. 



VI. ROBUST NON-LOCAL TWO-QUBIT GATE 

With high-fidelity local unitary gate and repeated 
QND measurement, the error probability for non-local 
coupling gates (e.g.. Fig. [s] and [4| is 



PC NOT (1 - F) + 2pL + 2s 



(11) 



which is dominated by the infidelity of the Bell pair 1 — F, 
since we assume ~ £m <^ l—F. In this section, we will 
show how to create high-fidelity Bell pairs between two 
registers with a reasonable time overhead. We will first 
briefiy review entanglement pumping [3TJ |35l [36] . Then 
we will quantitatively analyze the fidelity of the purified 
Bell pairs for our efficient two- level pumping scheme, and 
introduce the Markov chain model to calculate the fail- 
ure probability for entanglement pumping within a given 
number of attempts. Next we will quantify the perfor- 
mance of the high-fidelity Bell pair generation in terms of 
the total error probability (or average infidelity) and the 
time overhead, and discuss the trade-off between these 
two criteria. Finally, we will mention a non-post-selective 
pumping scheme which may further reduce the time over- 
head. 
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FIG. 6: (Color online) Two- level entanglement pumping be- 
tween registers and R-^ [circled by rounded rectangles in 
(a)]. (a,b) Generate and store one unpurified Bell pair, (c) 
Generate another unpurified Bell pair to pump (purify) the 
previously stored pair, (d) If the purification is successful, we 
obtain a purified Bell pair (level-1 pair) with higher fidelity; 
otherwise we discard the stored Bell pair and start the en- 
tire pumping process from the beginning, (e-f) The second 
level of entanglement pumping uses previously purified pairs 
to purify a stored Bell pair, to obtain a Bell pair with higher 
fidelity (level-2 pair). 



A. Entanglement pumping 

We now consider entanglement pumping [STJ |35l [36] 
with high-fidelity local unitary gate and robust measure- 
ment. During the entanglement pumping process, we 
first store one unpurified Bell pair [Fig.|6](a,b)], and then 
generate another unpurified Bell pair to purify the pre- 
viously stored pair [Fig. [6] (c,d)]. If the purification is 
successful, we will obtain a purified Bell pair with higher 
fidelity, which can be further purified by repeating the 
process in Fig. |6] (c,d) with new unpurified Bell pairs; 
otherwise we discard the stored Bell pair, and start the 
entire pumping process from the beginning. Sometimes, 
we may want to introduce a second level of entanglement 
pumping; that is to use previously purified pairs to purify 
a stored Bell pair [Fig. |6] (e-h)]. 



B. Fidelity of entanglement pumping 

We now analyze the performance of entanglement 
pumping for different errors of the unpurified Bell 
pairs. If the unpurified Bell pair is dominated by 
one type of error (e.g., dephasing error with den- 
sity matrix pdephasing = diag [F, 1 - F, 0, 0] in the 
Beh basis {|<l>+) , \^-) , |^+) , |^-)}, defined as \^^) = 
(|00) ± |11)) /a/2, = (|01) ± |10)) 7^2), we can skip 
the first level pumping. The unpurified pair then im- 
mediately becomes a level-1 pair and is purified with the 
circuit in Fig.jTjD. In Fig.jS^i we plot the fidelity curve (pu- 
rified fidelity v.s. number of successful pumping steps) 
for the one- level pumping process (i.e. 77,5 = 0), where a 
very high fidelity pair can be created after Up = 3 sue- 
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FIG. 7: (Color online) Bit-phase two- level entanglement 
pumping scheme to create high fidelity entangled pairs be- 
tween two registers and W . (a) Circuit for the first level 
pumping to purify bit-errors, corresponding to Fig. [g]:. (b) 
Circuit for the second level pumping to purify phase-errors, 
corresponding to Fig. |6^. The arrows indicate the time direc- 
tion for each register. Robust measurements are used here. 
If the two outcomes are the same, it is a successful attempt 
of pumping; otherwise we generate new pairs and restart the 
pumping operation from the beginning. 



cessful pumping steps. Note that we consider the full 
density matrix for all numerical calculation of entangle- 
ment fidelities [57], with the error models given in Eqs. 

If the unpurified Bell pair contains errors from 
both bit-flip and dephasing processes (e.g., depo- 
larizing error with density matrix ^depolarizing = 
diag [F, ^^^^ ^^^] -^^^^ basis), two- level en- 

tanglement pumping is needed. We introduce the follow- 
ing bit-phase two-level pumping scheme — the first level 
has rib steps of bit-error pumping using raw Bell pairs 
(Fig. 7^) to produce a bit-error-purified entangled pair, 
and the second level uses these bit-error-purified pairs for 



FIG. 8: (Color online) Entanglement fidelity F as a func- 
tion of the number of successful pumping steps, rib and Up 
are the number of pumping steps used to purify bit-errors 
and phase-errors, respectively. We assume fixed measure- 
ment and local two-qubit gate error rates sm — Pl — 10~^. 
(a) For bit-flip error pdephasing = diag [F, 1 — F, 0, 0] in the 
Bell basis { |$+) , |$") , |^+) , |^")}, with F = 0.95, one- 
level entanglement pumping is sufficient. High fidelity of 
Ffin = 99.98% can be achieved by rip = 3. (b) For de- 
polarizing error ^depolarizing = diag [F, ^^], with 
F = 0.95, two-level entanglement pumping is needed (see 
text for more details). The first level pumping only purifies 
the bit-error, but accumulates the phase-error at the same 
time, and therefore the (thin blue) fidelity curve for the ffist 
level pumping drops for rib > 1. The second level (thick red 
curve) uses the purified level- 1 pair (n^ = 1) to pump another 
stored pair. High fidelity of Ffin — 99.97% can be achieved 
by rip = 3. Note that we consider the full density matrix for 
all numerical calculation of entanglement fidelities [57], with 
the error models given in 



actor 1 — F, since measurement error does not cause infi- 
delity unless combined with other errors. In the limit of 



Up steps of phase-error pumping (Fig. [7)3). In Fig.jsjD we ideal operations {pl^Sm ^ 0), the infidelity e^'"^'''^^ 

vel (t 



-'£;,infid 



can 



plot the fidelity curves for the first level (thin blue curve) 
and the second level (thick red curve) of entanglement 
pumping. One-level pumping is insufficient to achieve 
high fidelity, but two-level pumping can achieve very high 
fidelity. With the parameters specified for Fig. [sJd, the 
maximum fidelity is achieved via the optimal choice of 
control parameters (n^,n*) = (1,3) for successful pump- 
ing steps of the first and second levels, respectively. 
For successful purification, the infidelity of the purified 
depends on both the control parameters 



(rib, Up) 

pair, e^; infid ' 



(77,5, np) and the imperfection parameters [F^pl^sm)- 
For depolarizing error, we find 

^(nb>i,np>i) _ 3 + 2np 4 + 2 (n5 + rip) 
^£;,infid ^ -A PL + (1 - j 



+ (np + 1) 



2(1-F) 



K + l)(l-F) 



(12) 



be arbitrarily close to zero, which is rigorously proved in 
Appendix |A] On the other hand, if we use the standard 
entanglement pumping scheme [35", "36] (that alternates 
purification of bit and phase errors within each pumping 
level), the reduced infidelity from two-level pumping is 
always larger than (1 — F)^ /9. Therefore, for very small 
Pl and ^m, the new pumping scheme is crucial to mini- 
mize the number of qubits per register. 

In Fis. 9 we show the contours of the infidelity infid^ 
as a funcTion of 715 and n^, where the contours are la- 
beled by values of log^^o ^^infid^ • parameters for the 
contour plots are F = 0.90 (left) and F = 0.95 (right); 
£m = PL = 10~^ (up) and Em = Pl = 10~^ (down). 
For optimal choice of (n^,n*), the minimal infidelity is 
limited by sm and pl- 

For dephasing error, one level pumping is sufficient (i.e. 
no bit-error purification, = 0). The infidelity is ap- 
proximately 



to leading order in pl and sm, for n^^rip > 1. The 
dependence on the initial infidelity 1 — F is exponen- 
tially suppressed at a cost of a linear increase of error 
from local operations pl and robust measurement Sm- 
Measurement-related errors are suppressed by the pref- 



.(0,np>l) 



(1 - F)""^^"^- 



-^PL + 2(1 -F)£M (13) 



by expanding to the leading order in pl and sm- 
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FIG. 9: (Color online) The contours of the infidelity £^ \^f^d 
as a function of rih and rip for depolarizing error. We use 

log^^Q £^ I'^f^^^ to label the contours. The other parameters 
are F = 0.90 (left) and F = 0.95 (right); sm = Pl = 10"^ 
(up) and sm — Pl — 10~^ (down). With optimal choice of 
(n^,n*), the minimal infidelity is comparable to the corre- 
sponding value of PL- 



VII. MARKOV CHAIN MODEL 

The overall success probability can be defined as the 
joint probability that all successive steps succeed. We 
use the model of finite-state Markov chain [58 to directly 
calculate the failure probability of (715, np)-two-level en- 
tanglement pumping using A/tot raw Bell pairs, denoted 

^(nf,,np) 



as s 



EJail 



Wot). 



A. Markov chain model for entanglement pumping 

We first use the Markov chain model to study the 
n-step one-level entanglement pumping. As shown in 
Fig. [lOj we use "0" to denote the initial state with no 
Bell pairs, "1" for the state with one stored unpurified 
pair, (j + 1) for the state with one purified pair surviving 
j steps of pumping, and " *" for the final state with the 
purified pair surviving n steps of pumping. Altogether 
there are n + 2 states. The (success) transition probabil- 
ity from state j to state j + 1 is g^, while the (failure) 
transition probability from state j to state is 1 — g^, 
for j = 0, 1, • • • , n. Here ^0 = 1 [corresponding to deter- 
ministic state transfer as shown in Fig. |6] (a,b)] and qj>i 




FIG. 10: (Color online) Markov chain model for one-level 
entanglement pumping. We use "0" to denote the initial state 
with no Bell pairs, " 1" for the state with one stored unpurified 
pair, (j + 1) for the state with one purified pair surviving j 
steps of pumping, and "*" for the final state with the purified 
pair surviving n steps of pumping. The (success) transition 
probability from state j to state j + 1 is qj , while the (failure) 
transition probability from state j to state is 1 — , for j = 
0, 1, • • • , n. Here ^0 = 1 and qj>i can be calculated according 
to the density matrix of the purified Bell pair surviving j steps 
of pumping 36 . The final state is self-trapped, and goes 
back to itself with unit probability. Each transition attempt 
consumes one unpurified Bell pair. 



can be calculated with the density matrix of the purified 
Bell pair surviving j — I steps of pumping The fi- 
nal state is self-trapped, and goes back to itself with unit 
probability, representing that once we have reached the 
desired final fidelity we no longer make any purification 
attempts and the system remains in this state with unit 
probability. Each transition attempt consumes one un- 
purified Bell pair. We would like to know the probability 
of reaching the final state " *" after A/tot attempts. More 
generally, we might also want to know the probability 
distribution over all n + 2 states. 

We use a (column) vector P with n + 2 elements to 
characterize the probability distribution among all n + 2 
states. From the t-th attempt to the {t + l)-th attempt, 
the probability vector evolves from Pt to Pt+i according 
to the rule 



Pt+i = M Pt, 
with the transition matrix 

/O 1-qi l-q2 ■■■ 1-qn 0\ 
1 

qi 

92 





(14) 



M 



(15) 



Since the initial probability vector is Pq = (1, 0, • • • , 0) , 
we can calculate the probability vector after A/tot at- 
tempts 



(16) 
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FIG. 11: (Color online) Markov chain model for two-level 
entanglement pumping. The required pumping steps are nt 
and Up for the two levels, respectively. We use "0, 0" to denote 
the initial state with no Bell pairs, "0, j + 1" for the state with 
one purified pair surviving j steps of pumping at the first level, 
"A; + 1, j + 1" for the state with one purified pair surviving 
k steps of pumping at the second level and one purified pair 
surviving j steps of pumping at the first level, "A; + 1, *" for 
the state with one purified pair surviving k steps of pumping 
at the second level and one purified pair surviving nt steps of 
pumping at the first level, and "*,0" for the final state with 
one purified pair surviving Up steps of pumping at the second 
level. For the first level pumping, the (success) transition 
probability from state "/c, j" to state "/c, j + 1" is qj, while the 
(failure) transition probability from state "/c, j" to state "/c, 0" 
is 1 — Qj, for j = 0, 1, • • • , rifo. For the second level pumping, 
the (success) transition probability from state "/c, *" to state 
"/c+l, 0" is Qfc, while the (failure) transition probability from 
state "/c, *" to state "0, 0" is 1 - Qk, for /c = 0, 1, • • • , Up. The 
final state is self-trapped, and goes back to itself with unit 
probability. 
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FIG. 12: (Color online) Failure probability s^;, fail as a function 
of A^tot- We assume a depolarizing error with F = 0.95, pl — 
£m = 10~^. We choose {ub^rip) = (2,3) for the lower curve, 
and (3,4) for the upper curve. For large A^tot, the failure 
probability Se, fan decreases exponentially with A^tot- 



The last element of pNtot success probability of 

reaching the final state "*,0" after A^tot attempts; the 
failure probability after A^tot attempts is then 



^(nb,np) 
-£;,fail 



(iVtc 



l-P{Ntot) 



(n5 + l)(np + l) + l • 



(18) 



In Fig. 

A^tot 
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we plot the failure probability s^E^aAi^ 







v.s. iVtot, lor control parameters {ub^rip) = (2,3) and 
(3,4). For A/tot sufficiently large, the failure probability 
decreases exponentially to zero. For any given param- 
eters, we can efficiently suppress the failure probability 
with some reasonably large A^tot- 



B. Total error probability &z average infidelity 



The probability vector PNtot describes the entire prob- 
ability distribution over all states of the Markov chain. 
The last element of P/v^^t is the success probability of 
reaching the final state "*" after A/tot attempts; the fail- 
ure probability after A^tot attempts is thus 



n+2 • 



(17) 



For two-level entanglement pumping, the state transi- 



tion diagram is shown in Fig. 11 and Up are the num- 
ber of pumping steps used to purify bit-errors and phase- 
errors, respectively. As detailed in Appendix [Bj we may 
use a (column) vector P with (n^ + 1) {rip + 1) + 1 ele- 
ments to characterize the probability distribution among 
all (77,5 + 1) {up + 1) + 1 states. From the t-th attempt 
to the {t -f- l)-th attempt, the probability vector evolves 
from P (t) to P (t + 1) according to the same rule as 
above [Eq. ([l4|)], but with the transition matrix M given 
in Eq. Q. 

Similar to one-level pumping, we can calculate the 
probability vector after A^tot attempts using Eq. (16). 
The probability vector pNtot describes the entire prob- 
ability distribution over all states of the Markov chain. 



We now introduce the total error probability (TEP) 
approximated by the sum of the failure probability and 
the infidelity of the purified Bell pair 



tot) 



^£;,fail 



(A^tot) 



„(nb,np) 
^£;,infid • 



(19) 



This is a very conservative estimate, since sometimes we 
do create some partially purified Bell pair though not the 
targeted purified Bell pair. And here we just say that the 
state has fidelity zero in these cases. 

To consider the possibility of using a partially purified 
Bell pair for output, we may introduce another useful 
quantity — the average infidelity (AIF) — for the output 
Bell pair from the robust entanglement generation, where 
we take into account these partially purified pairs. The 
average infidelity of the output pair is the weighted av- 
erage of the infidelity of the Markov chain 



= 1 



(20) 
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FIG. 13: (Color online) The optimized total error probability 
se [Eq. ( [21] )] (upper plots) and the optimized average infi- 
delity 5e [Eq. (22)] (lower plots) as a function A^tot- The er- 
ror probability for the local coupling gates is pi, = 10~^ (left 
plots) and pl — 10~^ (right plots). One-level pumping is used 
for dephasing error (thin blue curves); two- level pumping is 
used for depolarizing error (thick red curves). The other pa- 
rameters F — 0.95, Pi = pm — 5% are the same for all plots. 
Both £e and 5e saturate for large A^tot- In each plot, we also 
show the initial infidelity 1 — F (upper blue dashed lines) and 
the local error probability pL (lower violet dashed lines). 



Here the first term sums over all states of the Markov 
chain (except for the initial one), each of which has at 

least one partially purified pair with infidelity ^^j'nfid^ 
and probability ^ (^^tot)(nt,+i)n;+n;+2 5 ^^e last term 
comes from the situation that none of the partially pu- 
rified Bell pairs remain after the last attempt of the en- 
tanglement purification and we just use a classically cor- 
related pair with infidelity 1/2. Generally, the average 
infidelity is smaller than the total error probability. 

We may also optimize the choice of the control param- 
eters (n5,np) 



mm £ 



(iVtot) 



and 



(^tot) = ^^^^^^ ^E 



min S^:"'""^ {N,ot) 



(21) 



(22) 



In Fig. [13) we plot both the optimized total error prob- 
ability Ee and the optimized average infidelity Se as 
a function of A/tot- Both quantities asymptotically ap- 
proaches the same minimum value 



lim €E{Ntot) 

A^tot^oo 



lim Se (A^tot) = ^min- 

A/'tot^oo 



(23) 



Here the minimum value is simply the minimal infidelity 
of the entanglement purification 



Amin = min s 



{nh.rip) 
£;,infid 



(24) 



which is achieved by the control parameters 
(riij^np) = (nj,n*), for the imperfection parameters 
{PlA - F.em}' 

We remark that a faster and less resource intensive 
approach may be used if the unpurified Bell pair is dom- 
inated by dephasing error. Then, one-level pumping is 
sufficient (i.e. no bit-error purification, = 0). The 
optimized total error probability and average infidelity 
(thin blue curves) for this situation are plotted as a func- 
tion of iVtot in Fig. [HI 



C. Total time for robust entanglement generation 

The total time for robust entanglement generation Ie 
is proportional to the average number of raw Bell pairs 
generated (A^tot) 



iE-{Ntot) X {tE^tL^tM) 



(25) 



where t^; is the average generation time of the unpurified 
Bell pair. Note that the entanglement generation itself is 
a stochastic process. In principle, we may also include the 
stochastic nature of the entanglement generation by in- 
troducing a sub- level of Markov chain to characterize the 
stochastic entanglement generation. Since each entan- 
glement generation either succeeds or fails, the sub-level 
Markov chain only involves two states, which can be eas- 
ily incorporated into the Markov chain models discussed 
above. After incorporating the sub- level into the Markov 
chain, each transition corresponds to one attempt of en- 
tanglement generation, instead of one attempt of entan- 
glement purification that consumes one unpurified Bell 
pair previously. 

Nevertheless, the number of Bell pairs generated in a 
given period of time (i.e. A/tot) has a distribution. Since 

— 1/2 

the relative deviation of this distribution (~ ^tot ) 
fairly small for large A^tot (> 20), this only has a minor 
influence. Thus we replace (A^tot) by A^tot- 



D. Trade-ofF between gate quality and time 
overhead 



We now consider the balance between the "quality" of 
the robustly generated entangled pairs and the time over- 
head A^tot associated with the robust generation process. 
We may use either the optimized total error probability 
^E (Mot) or the optimized average infidelity Se (A^tot) 
to characterize the quality. Since both quantities ap- 
proaches the same asymptotic minimum Amin according 
to Eq. (23), there is only little improvement in the qual- 
ity of the robust entanglement generation once Ee (Mot) 
OT Se (Mot) is comparable to Amin (say 2Amin)- Thus, 
we find the value for A^tot by imposing the relations 



£e (Mot) = 2 A 



(26) 




FIG. 14: (Color online) Contours of the total error probability se (or average infidelity 6e) after purification (left), the total 
number of unpurified Bell pairs A^tot associated with se [Eq. (26)] (middle), and A^tot associated with 6e [Eq. (27)] (right). The 
contours are drawn with respect to the imperfection parameters pl (horizontal axis) and F (vertical axis). Two- level pumping 
(up) is used for depolarizing error, and one-level pumping (down) for dephasing error, pi = pM = 5% is assumed. 



or 



Se (iVtot) = 2A„ 



(27) 



First, we consider the total error probability Se (^tot)- 
The relation in Eq. (26) can be simplified, if we as- 



sume fixed control parameters (n^, Up) = (nj, n*) for the 
left hand side (rather than minimizing over all possible 
choices of {nij^np)). Combined with Eq. ([l9|, the failure 
probability should be comparable to the minimal infi- 
delity 



-£;,fail 



(A^tc 



(28) 



and the parameters (r 



, '6' '^pj 

F}, the above relation im- 
function a function of 



Since both the variable An 
depend on {pl,Pi,PmA 
plicitly determines A/tot 
{Pl^Pi^PmA- F}. 

In Fig. [m1 we plot the contours of Ee [Eq. ( [26] )] and 
A^tot [Eq. (|28p] with respect to the imperfection param- 
eters Pl and 1 — F, while assuming pi = pm = 5%. 
Actually the choice of pi and pm (< 10%) has negligi- 
ble effect on the contours, since they only modify Em 
marginally. For initial fidelity Fq > 0.95, the contours of 
Se are very close to vertical lines; that is £e is mostly 
limited by with an overhead factor (about 10) very 
insensitive to Fq. The contours of Ntot indicate that the 
entanglement pumping needs about tens or hundreds of 
raw Bell pairs to ensure a very high success probability. 



Similarly, we may also numerically obtain the value 
A/tot from Eq. (27). The contour plot of A^tot with re- 
spect to the imperfection parameters pl and 1 — F is 
also shown in Fig. [m] (c,f). We compare A^tots obtained 
from two estimates (total error probability [Eq. ( [26| )] and 
average infidelity [Eq. ([27[)]). As we expected, the A^tot 
obtained from total error probability is approximately 
1.2 ~ 2 time larger than the A/tot obtained from average 
infidelity, since the former is a more conservative estimate 
and requires more unpurified Bell pairs. Nevertheless, 
the difference is small and can be easily accounted by a 
prefactor of order unity. For clarity, in the rest of the 
paper we will use the A/tot estimated by using total error 
probability, and sometimes quote the values estimated by 
using average infidelity. 



E. Entanglement pumping with non-post-selective 
(NPS) scheme 

We now consider another entanglement pumping pro- 
tocol, proposed by Campbell [59 . The entanglement 
pumping scheme we have considered so far is post- 
selective (PS); that is we discard the Bell pair if one step 
of entanglement pumping is not successful. However, the 
Bell pair may still be highly entangled even if the entan- 
glement pumping failed at some intermediate step. The 
non-post-selective (NPS) entanglement pumping scheme 
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FIG. 15: (Color online) Markov chain model for one-level en- 
tanglement pumping with non-post-selective (NFS) scheme. 
The key difference from the previous Markov chain model 
with post-selective (FS) pumping scheme (see Fig. |lQ| ) is that 
here the transition for unsuccessful pumping reduces the chain 
label (score) by 1, while in the previous model the transition 
for unsuccessful pumping goes back to state "0" (restart of 
the entire pumping scheme). 




FIG. 16: (Color online) Markov chain model for two-level 
entanglement pumping with NFS scheme. We still use FS en- 
tanglement pumping scheme at the first level to have minimal 
accumulation of phase errors. Only at the second level, does 
the NFS scheme work more efficiently than the original FS 
scheme (see Fig. 11). 



[59] keeps track of the evolution of the density matrix 
of the Bell pair after each step of pumping. The NFS 
scheme avoids the inefficient restart (i.e., discarding in- 
termediately purified Bell pairs), and it may reduce the 
time overhead, especially when the unpurified Bell pairs 
have relatively low fidelity {F < 0.9). In Ref. [59 , the 
NFS pumping is discussed in the context of generating a 
graph state. 

We now describe how to use the NFS pumping scheme 
to generate purified Bell pairs. To simplify the discussion, 
we first assume that the errors from local measurements 
and operations are negligible. This assumption enables 
us to establish a connection between the Markov chain 
model and the NFS pumping scheme. 

Suppose the unpurified Bell pairs have only phase er- 
rors, then one level of entanglement pumping is sufficient. 
For this error model, one can show that a failed attempt 
produces an EFR pair with a density matrix identical to 
the one in the previous step [59 . One may introduce an 
accumulated score associated with entanglement pump- 
ing. The score increases by one unit for each attempt 
of successful pumping, and decrease by one unit for an 
attempt of unsuccessful pumping. The score for no Bell 
pair is 0, and for one unpurified Bell pair it is 1. The 
score exactly corresponds to the state label of the Markov 
chain (see Fig. [l5|. After each attempt of pumping, the 
score changes by ±1. If the score drops to (i.e. no 
Bell pair left), it gets back to 1 in the next attempt (i.e., 
creating a new unpurified Bell pair). The pumping pro- 
cedure continues, until the score reaches n + 1 (i.e., the 
final state "*" in the Markov chain). The key different 
from the previous Markov chain for post-selective pump- 
ing scheme (see Fig. 10) is that here the score decrease 
by 1 for unsuccessful pumping rather than restart from 
0. This modification increases the success probability of 
the robust entanglement generation. 



When the unpurified Bell pairs have both bit-fiip and 
phase errors (e.g., depolarizing error), we m ay use the bit- 
phase two- level pumping scheme (see Sec. VI B), which 



purifies the bit error at the first level and then the phase 
error at the second level. Since the phase error is not pu- 
rified at the first level, it accumulates after each attempt 
of pumping. Therefore, it is better to use FS entangle- 
ment pumping scheme at the first level to have minimal 
accumulation of phase errors. At the second level, the 
NFS scheme works more efficiently than the FS scheme. 
The Markov chain circuit for such mixed FS-NFS pump- 



ing schemes is shown in Fig. 16 



In practice, the error probability for the local oper- 
ations is always finite. Then our simple Markov chain 
model only provides an approximate description for the 
real process. The approximation comes from the fact that 
the score is now insufficient to specify the density matrix 
for intermediate Bell pairs, in the presence of local oper- 
ational errors. In order to obtain the density matrix for 
the intermediate state, we need to have the entire list of 
all previous pumping outcomes. Nevertheless, when the 
local operational errors are small compared to the infi- 
delity of the intermediate Bell pairs, the Markov chain 
model still provides an (optimistic) estimate for the total 
error probability and the average fidelity. 

We now compare the A/tot s associated with the FS and 
NFS schemes. The contours of the ratio between the two 
A^tots is plotted as a function of pl and F in Fig. [T7| As 
pointed out in Ref. [59], there is a significant improve- 
ment by using the NFS scheme (more than a factor of 3), 
for F < 0.9 and pl < 10""^. 
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FIG. 17: (Color online) The contours for the ratio between 
A^tot associated with the post-selective (PS) scheme and A^tot 
associated with the non-post-selective (NPS) scheme, as a 
function of and F. For both schemes, we use the same 
formula [Eq. (27)], but different Markov chain models (Fig. 10 
and |15[). The improvement from the NPS scheme becomes 
significant (more than a factor of 3), for F < 0.9 and pl < 



VIII. MAPPING TO DETERMINISTIC MODEL 

In this section, we will map our stochastic, hybrid, 
and distributed quantum computation scheme to a de- 
terministic computation model, which is characterized 
by two parameters — the clock cycle and the effective 
error probability. We will show that even when the un- 
derlying operations such as the entanglement generation 
are non-deterministic, our approach still maintains rea- 
sonable fast clock cycle time and sufficiently low effec- 
tive error probability. We will associate our discussion 
with achievable experimental parameters, consider the 
constraint set by the finite memory lifetime, and deter- 
mine the achievable performance of hybrid distributed 
quantum computation. 



A. Time and error in the theoretical model 

All the previous discussions can be summarized in 
terms of the clock cycle time 



tc — tE '^tL -\-tM^tE 
and the effective error probability 

7 EE ^2pL ^2£m, 



(29) 



(30) 



for a general coupling gate between two registers. 

We now provide an estimate of the clock cycle time 
based on realistic parameters. The time for optical ini- 
tialization/measurement is 



ti = tM 



\np 



M 



ln(l-7?)C' 



(31) 



with a photon collection/detection efficiency rj, vac- 
uum radiative lifetime r, and the cooperativity (Purcell) 
factor C for cavity-enhanced radiative decay [60l [61]. 
Eq. (31) is obtained from the estimate for the mea- 
surement error probability pm ^ ~ v) with 
^photon ^ ^m/ {'t/C). We assume that the entanglement 
is generated based on detection of two photons [3[ |32] , 
which takes time 



tE^{ti^r/C) /rj^ 



(32) 



Generally entanglement fidelity is higher for the two- 
photon schemes than one-photon schemes [31]. In addi- 
tion, some two-photon schemes have intrinsic purification 
against bit-flip errors [55 . The time for robust measure- 
ment is given in Eq. ([9]), and the total time for robust 
entanglement generation is given in Eq. ( [25] ). 

Combining Eqs. ([29]), (|3T|, (l32|, (l9| and (l25|, we ob- 



tain the clock cycle time (m units of the local operation 
time) as a function of other parameters 



tc 

tL 



f 



tLC 



(33) 



Meanwhile, we may obtain the relation m = 
^[Pl,Pi,Pm] by minimizing £m with Eq. ([8|, 
and find the relation A^tot = ^totiPL, F,sm] = 
Ntot [pl, 2Amin [pl^Pi^Pm]] using Eqs. ([24] and [26). 
Therefore, we have 



tc 

tL 



f 



tLC 



,PM,r],pL,F 



(34) 



The dimensionless parameter is the ratio between the 
times of emitting a single photon and performing a local 
unitary operation. For systems such as ion-traps and 
NV centers, this ratio is usually much less than unity 
(< 0.01). 

Similarly, we can obtain the effective error probability 
in terms of imperfection parameters 



1 = g[pL,PM,F], 



(35) 



by combining Eqs. ([8]), ([30]) and ([26|). 

In Fig. [18] we plot the clock cycle time tc and effec- 
tive error probabihty 7, for two-level pumping against 
depolarizing error. Assuming 77 = 0.2, we consider the 
two choices of parameters 1 — F = pi = pm = 5% (left) 
and 1% (right). For each case, we plot the contours of 
the normalized clock cycle time tc/tL as a function of 
Pl and and the effective error probability 7 as a 

function of pl- The clock cycle time can be reduced by 
having a fast radiative decay rate r/C, which can be fa- 
cilitated by having a large cooperativity factor C. The 
reduction of the clock cycle time stops once this ratio is 
below certain value, approximately 0.003 (left) and 0.001 
(right), where local gate operation becomes the dominant 
time consuming step. Similarly, we plot the clock cycle 
time tc and effective error probability 7, for one-level 
pumping against dephasing error in Fig. [19] 



14 




FIG. 18: (Color online) Plots of clock cycle time tc and ef- 
fective error probability 7, for two- level pumping against de- 
polarizing errors. Upper plots: contours of the normalized 
clock cycle time tc/th (upper) as a function oi pl and y~c 
(the normalized effective radiative lifetime). Lower plots: the 
effective error probability 7 as a function oi pl- We assume 
1 - F = PI = pM = 5% (left) and 1% (right), and ry = 0.2. 



FIG. 19: (Color online) Plots of clock cycle time tc and ef- 
fective error probability 7, for one-level pumping against de- 
phasing error. Upper plots: contours of the normalized clock 
cycle time tc/th (upper) as a function oi pl and (the 
normalized effective radiative lifetime). Lower plots: the ef- 
fective error probability 7 as a function oi pL- We assume 
I - F = PI = pM = ^% (left) and 1% (right). And r] = 0.2. 



In the limit of negligible radiative decay time, we ob- 
tain the lower bound for the normalized clock cycle time 

i^w. + /+ / ^tot for m = . . 

where for m > 1 there is a time overhead 2m + 2 associ- 
ated with local operation and robust measurement; while 
there is no such overhead for m = 0. 



memory error probability for each clock cycle is approxi- 
mately tc/tmem^ which should be small (say 10~^) in or- 
der to achieve fault-tolerant quantum computation. This 
is indeed the case for both trapped ions (where t^em ^ 10 
s has been demonstrated [U, ^47 ) , as well as for proximal 
nuclear spins of NV centers (where tmem approaching a 
second can be inferred [H]). So far, we have justified the 
feasibility of the hybrid distributed quantum computa- 
tion scheme. In the next subsection, we will provide a 
criterion for hybrid distributed quantum computation. 



B. Estimated numbers for experimental setups 

Suppose the parameters are (ti,,r, ?7,C) = 
(0.1 /is, 10 ns, 0.2, 10) iSl ESI and 
{I-F.pi.pm^Pl.sm) = (5%,5%,5%,10-^8x 10-4) 
for our quantum registers (based on ion-traps or NV 
centers). For depolarizing errors, two- level pumping 
can achieve (tc,7) = (200 /is, 2.7 x 10"^) . For some 
entanglement generation schemes jSj E2l principle 
only dephasing error exists, because they have intrinsic 
purification against bit-flip errors. If all bit-flip errors 
are suppressed, then one-level pumping is sufficient and 
(tc,7) = (42 /is, 2.2 X 10-^)[74 . In Table |lj we have 
listed (tc7) for parameters 1 — F = pm = Pi = 5% or 
1%, and PL = 10-^, 10-"^, 10-^ or 10"^. As expected, 
we find that tc gets longer, if the fidelity F is lower 
and/or the error probability {pm or pi) is higher; tc 
significantly reduces if the error for the unpurified Bell 
pairs changes from depolarizing error to dephasing error. 

We remark that tc should be much shorter than the 
memory time of the storage qubit, tmem- Because the 



C. Constraints from finite memory life time 

Above we have mostly ignored the effect of finite mem- 
ory time, and with the various sequences of purification 
of imperfections the final fidelity of the operations have 
then been limited only be the local operation. All of 
these purifications, however, increase the time of the op- 
erations and eventually the system may become limited 
by the finite life time of the memory. In this subsection 
we shall evaluate this constraint set by the finite memory 
lifetime. 

To simplify the discussion we assume that we have a 
very short radiative lifetime r or that we are able to 
achieve a very large Purcell factor so that r/C becomes 
negligible. All the time scales are then proportional to 
the local gate time tL- With a finite memory time, i.e., 
some fixed tmem/tL^ there is a limit to have many oper- 
ations we can do before we are limited by the memory 
error. To get an estimate for this limit we assume that 
the ideal number of operations is roughly given by the 



15 









Depolarizing 








Dephasing 










F 


= 0.95 


F 


= 0.99 




F 


- 0.95 


F 


= 0.99 










7 




7 






7 




7 




VL 


= 10"^ 


65 


1.9 X 10"^ 


19 


9.1 X 10" 


-3 


20 


1.7 X 10"^ 


8 


8.7 X 10" 


-3 


PL 


= 10-^ 


200 


2.7 X 10"^ 


49 


1.2 X 10" 


-3 


42 


2.2 X 10"^ 


17 


9.9 X 10" 


-4 


PL 


= 10-^ 


387 


3.5 X 10-^ 


65 


1.3 X 10" 


-4 


80 


2.8 X 10"^ 


22 


1.2 X 10" 


-4 


PL 


= 10-^ 


997 


4.5 X 10"^ 


162 


1.7 X 10" 


-5 


140 


3.4 X 10"^ 


39 


1.4 X 10" 


-5 



TABLE I: We list the values of tc and 7 as a function of (rows) and F (columns) for depolarizing and dephasing errors of 
the unpurified Bell pairs. We also assume pm — pi — 1 — F ^ and (tL^r, 77, C) = (0.1 //s, 10 ns, 0.2, 10). Note that tc estimated 
by using average infidelity is approximately 1.3 ~ 1.6 times less than the numbers listed here. 



point, where the memory error probability is the same as 
the effective error probability for the non-local coupling 
gate: 



^ 7- 



Then according to Eq. ([36|), we have 



tc 



tL 



tc tL 



(37) 



= 7-^(2m + 2-(^^,o)A^tot, (38) 



where the variables {7, m, A^tot} are all determined by the 
imperfection parameters {1 — F^pm^Pi^Pl}- We further 
reduce the imperfection parameters by assuming 1 — F = 
Pm = Pi^ and get the contour plot of tmem/tL in terms of 
the imperfection parameter pl and 1 — F in Fig. [20] In 
the plot, we consider both the situation of depolarizing 
or dephasing error during entanglement generation. 

For given tmem/tL, we may use Fig. 20 to find the 
valid region in the parameter space of pl and 1 — and 
then identify the achievable effective error probability 7. 
For example, with ion-trap systems it may be possible 
to achieve tmem/tL ~ 10^ [43, 44| [62], and the region 
left of the shaded contour line ( log ^ q ^mem/^L = 8) can 
then be accessed, which enables us to obtain a wide range 
effective error probability 7 depending on the practical 
values of PL and 1 — F. For NV centers, it should be 
feasible to achieve tmem/tL ^ 10^ by having tmem ^ 10 
s and tL ^ 10~^ s [12 ; the region on the left side of 
the shaded contour line {log^Q tmem/tL = 7) still covers 
a large portion of the parameter space. For a given ex- 
periment situation with a finite memory time as well as 



other imperfections, we can thus use Fig.[20]to determine 
the achievable performance of hybrid distributed quan- 
tum computation. 



IX. APPROACHES TO FAULT TOLERANCE 

The entanglement based approach discussed in this pa- 
per provides a method to make gates between any quan- 
tum registers and this can be used to implement arbi- 
trary quantum circuits, once the errors in the gates are 
sufficiently small. The errors can be further suppressed 
by using quantum error correction. For example, as 



shown in Table [l| {pljF) = (lO ^, 0.95) can achieve 7 ~ 
2.7 X IQ-^, well below the 1% threshold for fault toler- 
ant computation based on approaches such as the C^/Cq 
code [65| or 2D toric codes ^66]; {pl^F) = (10-^,0.99) 
can achieve 7 1.7 x 10~^, which allows efficient codes 
such as the BCH [[127,43,13]] code to be used without 
concatenation. 

Following Ref. p7] we estimate 20 registers per logical 
qubit to be necessary for a calculation involving = 10^ 
logical qubits and Q = 10^ logical operations, assuming 
the memory failure rate and effective error probability are 
tc/tmem ^ 7 ^ 1-7 X 10"^ (e.g., achieved by tmem ^ 10 
s, tc ^ 162 /is). (This estimate is based on Fig. 10b 
of Ref. [63.) Assuming that error correction is applied 
after each logical operation, and that logical operation 
and following recovery take approximately 2 — 16 clock 
cycles depending on the type of operation and the coding 
scheme (see section II. A of Ref. [67 ), the total running 
time of this computation would then be approximately 
400 - 3000 s. 

We remark that one important property of distributed 
quantum computation is that the measurement time is 
relatively fast compared with the non-local coupling gate, 
because the measurement does not rely on the time- 
consuming processes of entanglement generation and pu- 
rification while the non-local coupling gate does. This 
property is different from the conventional model of 
quantum computation, where the measurement is usu- 
ally a slow process that induces extra overhead in both 
time and physical resources [67^. Thus, instead of recon- 
ciling slow measurements [68] , it might also be interesting 
to study possible improvement using fast measurements 
for fault-tolerant quantum computation. 

The above estimates have been performed assuming 
that our hybrid register based approach is mapped di- 
rectly to the standard circuit model. In some situations 
this may, however, not be the most advantageous way 
to proceed, since the register architecture may allow for 
more efficient performance of certain tasks. As a par- 
ticular example, we now briefly discuss an alternative 
new approach to fault-tolerant preparation of GHZ states 
(e.g., |00 • • • 0) + |ll • • • 1)), which are a critical component 
both for syndrome extraction and construction of univer- 
sal gates in quantum error correcting codes [65l |69] . This 
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FIG. 20: (Color online) Contours for tmem/tL (shaded contours) and 7 (red contours), as a function of pi, and 1 — F. The 
label for the contour values are in logarithmic scale with base 10. We consider both situations of (a) depolarizing error and (b) 
dephasing error. The other parameters are pM — pi — 1 — F and t /C — 0. 
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- — 



FIG. 21: (Color online) Circuits for fault-tolerant prepara- 
tion of GHZ state, (a) Conventional circuit uses 8 quantum 
registers to generate a 4-qubit GHZ state. The state |+) is 
(|0) + |1)) /a/2, (b) New circuit requires only 4 registers, by 
using partial Bell measurement (PBM) between registers, de- 
tailed in (c). 



new approach relies upon the observation that the EPR 
pairs from entanglement generation can be used for de- 
terministic partial Bell measurement (PBM), which is 
achieved by applying local coupling gates and projective 
measurements as shown in Fig. [21) 3. After the PBM, 
the two storage qubits are projected to the subspace 
spanned by the Bell states \^^) = (|00) ± |11)) /V^ if 
the measurement outcomes are the same, or they are 
projected to the subspace spanned by the Bell states 
1^^) = (|01) ± |10)) /a/2 if the measurement outcomes 
are different. For the latter case, we may further flip one 
of the storage qubits, so that they are projected to the 



subspace spanned by |$^). By using PBMs, we can per- 
form fault-tolerant preparation of GHZ state efficiently 
(up to single qubit rotations) as detailed below. 

Fault-tolerant state preparation requires that the prob- 
ability to have errors in more than one qubit in the pre- 
pared state is O (p^), with the error probability for each 
input qubit or quantum gate being O (p); that is multi- 
ple errors only occur at the second or higher orders. The 
regular circuit to prepare a four-qubit GHZ state fault- 
tolerantly ^5] is shown in Fig. 21 i. If this circuit should 



be implemented with quantum registers, the CNOT gates 
in Fig.f2l|i should be created by using the circuit detailed 
in Fig. [3] and eight quantum registers would be required. 

By using PBMs, however, only four quantum regis- 
ters are needed in order to generate GHZ states fault- 
tolerantly as shown in Fig. [STJd. The fault-tolerance 
comes from the last (redundant) PBM between the sec- 
ond and fourth register (Fig. [21p), which detects bit- 
errors from earlier PBMs. The advantage of the PBMs 
is that it propagates neither bit- nor phase-errors. The 
circuit of Fig. 21 z indicates that the only way to propa- 
gate error from one input to another (say. Si to S2) is via 
some initial error in the entangled pair between Ci and 
C2. However, for Bell states |^^), we have the following 
identities 



(39) 
(40) 



Similar identities also exist for Bells states |^^). Suppose 
Si has an error, because of the above identities, we can 
always treat the imperfection of the Bell pair |$+)^^ 
as an error in Ci (the qubit from the same register as 
Si). Therefore, only the first register has errors and they 
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FIG. 22: Circuit for fault-tolerant preparation of the 
GHZ state with only 8 quantum registers. 




-qubit 



never propagate to the second one. (PBM may induce 
errors to unconnected but entangled qubits.) 

The present scheme may be expanded to larger num- 
bers of qubits and generally we may fault-tolerantly pre- 
pare 2^-qubit GHZ state with only 2^ quantum registers, 
by recursively using Fig. [21) 3 with the two dashed boxes 
replaced by two 2^~^-qubit GHZ states. The circuit for 
fault-tolerant preparation of the 8-qubit GHZ state is 
shown in Fig. 22 Note that we can perform PBMs act- 



ing on different registers in parallel. Suppose each PBM 
takes one clock cycle, the preparation time is 2 (clock cy- 
cles) for a 4-qubit GHZ state shown in Fig. 



21 



3. The two 



PBMs in the orange boxes are performed in the first clock 
cycle, and the rest for the second clock cycle. Generally, 
for a 2^-qubit GHZ state with n > 3 (see discussion in 
Appendix [C]) , the preparation time is only 3 (clock cy- 
cles), and the error probability for each register is only 
approximately 3p/2. Therefore, the PBM-based scheme 
for fault-tolerant preparation of the GHZ state is efficient 
in both time and physical-resources. 
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FIG. 23: (Color online) The contours of the infidelity £^Ei^d^ 

as a function of m and rip. We use logio ^L^infid^ label 
the contours. We assume a depolarizing error with initial 
fidelity F = 0.95, and sm — Pl — 0. The final infidelity 
can be arbitrarily small for sufficiently large rih and rip. This 
indicates that the bit-phase two-level entanglement pumping 
scheme can create pairs with arbitrarily high fidelity. 



ble to further facilitate fault-tolerant quantum compu- 
tation with systematic optimization using dynamic pro- 
gramming [ 7Qj . 
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X. CONCLUSION 

In conclusion, we have proposed an efficient register- 
based, hybrid quantum computation scheme. Our 
scheme requires only five qubits (or less) per register, and 
it is robust against various kinds of imperfections, includ- 
ing imperfect initialization/measurement and low fidelity 
entanglement generation. We presented a Markov chain 
model to analyze the time overhead associated with the 
robust operations of measurement and entanglement gen- 
eration. We found reasonable time overhead and consid- 
ered practical implementation of quantum registers with 
ion traps or NV centers. We also provided an example 
using partial Bell measurement to prepare GHZ states for 
fault-tolerant quantum computation. It might be possi- 



APPENDIX A: BIT-PHASE TWO-LEVEL 
PUMPING SCHEME 

In this Appendix we show that our bit-phase two-level 
entanglement pumping scheme can create pairs with fi- 
delity arbitrarily close to unity, if we have perfect lo- 
cal operations. Numerical indication of this is shown in 
Fig. [23j and in the following we will provide a rigorous 
proof to this claim. 

We assume that the initial state is a mixed state that 
has only diagonal terms in the Bell basis 

p = a (<l>+|+6 \^-) (<!>" l+c |^+) 1^-) 

(Al) 
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where = (|00) ± |11)) /a/2, I^-*) 

(|01) ± |10)) /V2, and the coefficients are non-negative 
and sum to unity. (This assumption is only made to 
simplify the presentation. For a general density matrix, 
only the diagonal elements given in Eq. (Al) are im- 
portant [34 .) After the purification the density matrix 
retains this form but with new coefficients. Therefore, 
we only need four coefficients for each state using Bell 
basis, denoted as the fidelity vector F = (a^b^c^d). 
We use the lower index to keep track of the pumping 
steps, so that the fidelity vector for the unpurified 
state is Fq = (ao, 60, cq, (io), and the vector for the 
purified state after n steps of entanglement pumping is 

Suppose we use the Fq state to pump the state F^ 
against bit-errors, the success probability is 

Pn+l = {CLO + ^0) {cLn + K) + (cq + do) (Cn + dn) , (A2) 

and the fidelity vector is 



1. First level pumping 

For the ffist level of pumping against bit-errors, we 
may parameterize the fidelity vector 



F^n — is^ni ^n) 



(A3) 



2 ■ 2 

(AlO) 

in terms of two variables 5n and rjn^ and obtain some 
bounds for these variables. Since we are pumping against 
bit-errors, r]n decreases with n 

Vn+i <r]n<-- - <m = a/2< 1/6. (All) 

The success probability for the (n + l)th step of pumping 
is 

Pn+i = (1 - a) (1 - 2r]n) + 2a7?n (A12) 
< 1 - a - 2ar]n < 1 - a, (A13) 

and on the other hand we have 

Pn+i > (1 - a) (1 - 2r]n) >^(l-a) (A14) 
where the second inequality follows from rjn < 1/6 [74] . 



{aoan + bobn, a^b^ + ^oa^, cqc^ + d^dn, c^dn + d^Cn) , The third element of the fidelity vector is 



(A4) 



Pn+1 



for perfect local operations (measurement and CNOT 
gate). 

Similarly, for pumping against phase-errors, the suc- 
cess probability is 

Pn+l = i^O + Co) {an + Cn) + (^0 + ^^o) {K + dn) , (A5) 

and the fidelity vector is 

(A6) 



Since Pn+i > I (1 — o^), we have 



< 



3 (1 - a) 

3X n+l 

4 



Vn < 



< 



Pn+1 



3a 



(A15) 



2(1 - a) 



n+l 



CO (A16) 



Co 



Pn+1 



{aoan + CoCn, bobn + dodn, aoCn + Cottn, bodn + dobn) . 

(A7) 



which indicates that r]n approaches zeros exponentially 
with respect to n. When n is large enough, r]n is negligi- 
ble and Pn ^ I — — O ijin)' 

Similarly, we obtain the recursive relation for 6n- 



In general, any state can be turned into a so-called 
Werner state with the same fidelity, and as a worst case 
scenario we shall assume the unpurified Bell state to be 
a Werner state: 



1 1 — a — 2ar]n ^ 1 — 2a 

1 



Fq — ( Fq 



1 — i^o 1 — i^o 1 — i^o 



Thus 



(A8) 



with Fq > 1/2 to ensure that it contains distillable entan- 
glement. For convenience of later discussion, we rewrite 
Fq as 



> 



1 -2a 



1 -2a 
1 — a 



(A17) 
(A18) 
(A19) 



1 



Sn> - > 



I -2a 



1 — a 



n+l 



(5o (A20) 



l-2a\''+^ l-3a 



Fo = (1 - 3a/2, a/2, a/2, a/2) 



(A9) 



and 5q 



where a = | (1 - Fq) < |. 



1 — a 



^ . We have therefore 

^ < p-4a \" 

Sn V 3a / 1 — 3a 



a 



(A21) 



19 



2. Second level pumping 



If i > > ^0 > > Vo^ we have 



After rib steps of first level pumping, we have the fi- 
delity vector 



(A22) 



/ \ n-t-i 

where Sq = 5n^ > (^T^j and r/g = r/n^ < 



, n+1 



Co. 



The fidelity vector after Up steps of pumping against 
phase errors is 

Fl,= {]^+^'n,\-5'n-ri'n- dn, v'n, dn^ , (A23) 

where 6'„<l/2 and t]'^ < 1/2. 

The success probabihty for the (n + l)th step is 



Pn+l 



1 



1 



0-Vo 



= 2+2('5o+^o)('5;+r?;) 
1 

>2 



(A24) 



We now consider the elements of F^. On one hand, the 
erroneous admixture of |^~) described by {dn} keep de- 
creasing with n, since |^~) errors are also purified during 
the second level pumping. On the other hand, the erro- 
neous admixture of |^^) described by {rj!^} may increase 
with n, but it is upper-bounded by the following relation: 



Vn+i — 



aoCn + Co an 



< 2 



Pn-\-l 

<(l + 2^^r+%^. 



and one can also show the lower bound for 77^+ 1 by in- 
duction: 



However, S!^ approaches 1/2 much faster: 



(A26) 



+ — 

Pn+l 



(A27) 



2 ■ l + 46'oS'^+4{S'oV'n+5M+vM 



+ 



^Wn + S'nV'o + vWn < ^Wn + V'o < S'oC, (A28) 

where we in the second inequality introduced a number 
C such that SqT]^^ -\- tjq < Sq( . In the next subsection, we 
will show that we may choose (" small (i.e., (" = 2s) such 
that 



a: 



n+l 



> 



6'„ + 5',{l-2C) 
1 + 45^,5; (1 + C) 
a; + (1 - 20 



i + 4a;,(i-2C)a; 

>A:, + AU1-2C)-2AU1-2C)A:, 



(A29) 



(A30) 



where A^, = ^x^(5; « (1 + 3C/2) S'^ for n = 0,1, ■ ■ ■ , 
and the third equality uses {^^y > x + y — 2xy for 
< x,y < 1/2. Finally, we have 

1/2 - A;+i < (1 - 2X'o (1 - 20) (1/2 - a;) 

<(l-2AUl-20r^'(l/2-A;,) (A31) 

So far, we have introduced inequalities to bound ele- 
ments of fidelity vectors at different stages of pumping. 
In the next subsection, we will use these constraints to 
show that we are able to achieve fidelity arbitrarily close 
to unity by carefully choosing the numbers of pump- 
ing steps {rib, rip) for the bit-phase two- level pumping 
scheme. 

3. s — N argument 

For > and a < 2/7, we may choose 



rib ^ naax • 



In £ — In -. 



31ne-ln - 



In . 



In 



2(l-«) 



(A25) such that 



Vo < 



3a 



< 



2{l-a)J 2 



3a 



and 



2-Aa I -3a 



1 -2a\ ' I -3a 



< e 



I- a 



Then we choose for n„ 



21n£ 



ln(l-25i) >'^P> ln(l-2AJ,(l-2C))- 



In £ 



(A32) 



(A33) 



(A34) 



(A35) 



(A36) 
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Such n„ always exists, for ^ ^ > 
, ^ -y. Thus, we have 

-25^,^(1+0(1-20 ln(l-2AJ,(l-20) 

7?;^ < (1 + 2J^)"- < (1 - 2J^)-"^ vo < e-% < ^0 

(A37) 

and 

1/2 - < (1 - 26^, (1 - 20)"^ (1/2 - A[,) < £/2 

(A38) 

We now verify the two required relations which are 
underlined in the previous discussion. The relation 
^ > > ^0 ^ Vn > Vo is satisfied for all < n < n^. 
And the relation SqT]'^^ + < ^oC is also implied, since 

Finally, the achievable fidelity for the above choice of 
ni) and Up is 

1 - F„^,„, = 1/2 - 5'^^ = 1/2 - (j^) 

< 1/2 - (\^) (1/2 - £) < 4e. (A39) 

The bit-phase two-level pumping thus allows us to ap- 
proach F = 1 with arbitrary good precision. 

APPENDIX B: MARKOV CHAIN MODEL FOR 
TWO-LEVEL PUMPING 

Here we present the Markov chain model for two-level 
entanglement pumping. 

The state transition diagram for two-level entangle- 
ment pumping is shown in Fig. [TT] We assume that the 



required pumping steps are 77.5 and Up for the two lev- 
els, respectively. Since two entangled pairs are stored, 
we need two labels to track the intermediate state for 
two-level entanglement pumping. 

Here we use "0,0" to denote the initial state with no 
Bell pairs, "0, j + 1" for the state with one purified pair 
surviving j steps of pumping at the first level, "/c+1, j+1" 
for the state with one purified pair surviving k steps of 
pumping at the second level and one purified pair sur- 
viving j steps of pumping at the first level, "/c + 1,*" 
for the state with one purified pair surviving k steps of 
pumping at the second level and one purified pair surviv- 
ing rib steps of pumping at the first level, and "*,0" for 
the final state with one purified pair surviving Up steps 
of pumping at the second level. 

For the first level pumping, the (success) transition 
probability from state "/c, j" to state "/c, j + 1" is qj, 
while the (failure) transition probability from state "/c, j" 
to state "A:,0" is 1 — qj, for j = 0, 1, • • • ,715. For the 
second level pumping, the (success) transition probabil- 
ity from state "/c, *" to state "/c + 1,0" is Qfe, while the 
(failure) transition probability from state "/c,*" to state 
"0, 0" is 1 — Q/e, for /c = 0, 1, • • • , n^. Here the transition 
probabilities {qj} and {Qk} can be calculated accord- 
ing to the density matrices of the intermediate purified 
Bell pairs [36 . The final state is self-trapped, and goes 
back to itself with unit probability. Altogether there are 
(77,5 + 2) {up + 1) + 1 states. 

In order to fulfill the requirement that each transition 
attempt consumes one unpurified Bell pair, we need to 
contract the states of "/c, *" and "/c + 1, 0" into one state, 
since this transition does not consume any unpurified Bell 
pair. After the contraction, there are (n^ + 1) {up + 1) + 1 
states remaining. 



Therefore, we may use a (column) vector P with (n^ + 1) {up + 1) + 1 elements to characterize the probability 
distribution among all (n^ + 1) {rip + 1) + 1 states. From the t-th attempt to the {t + l)th attempt, the probability 
vector evolves from P (t) to P {t -\- 1) according to the following rule 



P(t + l) = MP(t), 



and the transition matrix is 



M = 



/Mi+^n,(l-go)Ni gn,(l-Qi)Ni ••• gn,(l-Q2)Ni 0\ 
^n.QoNi Ml 



\ 







Ml 



with sub-matrices: 



Ml 



/ 1-qi 
1 

Qi 

\ 







(nb + l)x(nb + l) 



(Bl) 



(B2) 



(B3) 
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and 



/o •• 



Ni = 



1 \ 




(B4) 



\o J 



(n{, + l)x(nf, + l) 



APPENDIX C: FAULT-TOLERANT 
PREPARATION OF 2^-QUBIT GHZ STATE 

We consider fault-tolerant preparation of a 2^-qubit 
GHZ state with n > 3. We label 2^ registers by 
0, 1, • • • ,2^ — 1. The GHZ state can be prepared in just 
3 (clock cycles): in the first clock cycle, we perform PBM 
for register pairs (2m, 2m + 1) with integer m; in the sec- 
ond clock cycle, we perform PBM for pairs (2m, 2m + 2) 
and (2m -h 1, 2m + 3); in the last clock cycle, we per- 
form PBM for pairs (8m + 3, 8m + 5), (8m + 4, 8m -j- 6), 
(16m' -h 7, 16m' + 9), (16m' + 8, 16m' + 10), 

(2^-1 + 2^-2 - 1, 2^-1 + 2^-2 + 1) , and 
(2^-1 + 2^-2, 2^-1 -h 2^-2 -h 2) in parallel. 

In order to prepare the specific GHZ state |00 • • • 0) -h 



|11 • • • 1), we still need to perform bit-flip operations for 
individual registers, which are determined by the mea- 
surement outcomes for all PBMs. Suppose the error 
probability for each PBM is p. Since the redundancy 
checks of the PBMs impose consistency requirements for 
measurement outcomes (error detection) , the probability 
for undetected errors in measurement outcomes has been 
suppressed to O (p^) for each PBM. To the leading order 
of j9, we only need to consider the phase errors from PBM 
that are not detected by the redundancy check. Thus, 
the total error probability scales as 2^~^3p, and the er- 
ror probability for each register is only approximately 
3p/2. Therefore, we have demonstrated a PBM-based 
scheme to prepare the GHZ state fault-tolerantly, which 
is efficient in both time and physical- resources. 
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